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Cosmic rays with energies beyond the Greisen-Zatsepin-Kuzmin 'cutoff' 
at ~ 4 x 10 10 GeV pose a conundrum, the solution of which requires either 
drastic revision of our astrophysical understanding, or new physics beyond 
the Standard Model. Nucleons of such energies must originate within the lo- 
cal supercluster in order to avoid excessive energy losses through photopion 
production on the cosmic microwave background. However they do not point 
back towards possible nearby sources, e.g. the active galaxy Cen A or M87 in 
the Virgo cluster, so such an astrophysical origin requires intergalactic mag- 
netic fields to be a hundred times stronger than previously believed, in order 
to isotropise their arrival directions. Alternatively the primaries may be high 
energy neutrinos, say from distant gamma-ray bursts, which annihilate on the 
local relic background neutrinos to create "Z-bursts". A related possibility 
is that the primary neutinos may initiate the observed air showers directly if 
their interaction cross-sections are boosted to hadronic strength through non- 
perturbative physics such as TeV-scale quantum gravity. Or the primaries 
may instead be new strongly interacting neutral particles with a longer mean 
free path than nucleons, coming perhaps from distant BL-Lac objects or FR-II 
radio galaxies. Yet another possibility is that Lorentz invariance is violated 
at high energies thus suppressing the energy loss processes altogether. The 
idea that has perhaps been studied in most detail is that such cosmic rays 
originate from the decays of massive relic particles ( "wimpzillas" ) clustered 
as dark matter in the galactic halo. All these hypotheses will soon be criti- 
cally tested by the Pierre Auger Observatory, presently under construction in 
Argentina, and by proposed satellite experiments such as EUSO. 
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1 Introduction 



In 1962, Linsley JIJ detected a cosmic ray air shower at the pioneering Volcano Ranch 
array with an estimated energy exceeding 10 11 GeV. The significance of this event became 
clear only some years later, following the serendipitous discovery of the universal 2.7 K 
cosmic microwave background (CMB), when Greisen and, independently, Zatsepin and 
Kuzmin || noted that cosmic ray nucleons with energies exceeding Eqzk ~ 4 x 10 10 GeV 
would suffer severe energy losses through scattering on the CMB photons due to excitation 
of the A (1232) resonance with a cross-section of ~ 400 fib. Consequently their typical 
range should decrease rapidly above this energy, leading to a 'GZK cutoff' in the energy 
spectrum if the sources are cosmo logically distant. Conversely the absence of such a cutoff 
in the spectrum (which is approximated by a power-law dN/dE cx E~ 3 at lower energies) 
would imply that the sources are nearby, within the local supercluster.f] 

Plausible astrophysical sources for such ultra-high energy cosmic rays (UHECRs) were 
discussed by Hillas || who noted that simply to contain charged particles long enough 
for them to be accelerated to such high energies required sufficiently large magnetic fields 
and/or spatial volumes such as could plausibly be achieved only in a few specific sites such 
as the extended lobes of giant radio galaxies or active galactic nuclei. Moreover particles of 
such energies should be undeflected by the weak intergalactic magnetic fields, so should 
point back to their sources. Thus there was the exciting prospect that further such 
observations would clarify the long-standing question of the origin of cosmic rays. Indeed 
UHECRs continued to be observed with other detectors, notably at Haverah Park || and 
Yakutsk JTO] . However, perhaps because of the theoretical prejudice that the sources of 



UHECRs must be cosmologically distant, as well as concerns regarding uncertainties in 
the energy measurements, these observations appear not to have been taken seriously JTT 



As late as 1993 it was stated that the GZK cutoff is indeed present in the data, confirming 



the expectation that the sources of UHECRs are distant FR-II radio galaxies [12[ . 

The breakthrough came the same year with the detection of an event of energy 
(3.2 ± 0.9) x 10 11 GeV by the innovative Fly's Eye detector in Utah, which measured 
the fluorescence from the excited N 2 molecules in the atmospheric shower and was thus 



able to reliably determine the energy [|13f| . Moreover the shower was seen to have its 
maximum at an atmospheric depth of 815^35 gem -2 , consistent with a proton primary, 
but significantly less than the expected value of 1075 gem -2 for a photon primary fl4||f1 
Subsequently the Akeno Giant Air Shower Array (AGASA) has detected a large number 
of air showers with energies exceeding E GZK [[L7|] . 

* Recent detailed calculations Q indicate that the typical range of a nucleon drops below ^100 Mpc 
above 10 11 GeV. For heavy nuclei the dominant energy loss process is photodissociation; above 10 11 GeV 
this happens through scattering on the CMB and the typical range is even smaller than for nucleons Q . 
For photon primaries the dominant opacity above Eqzk is due to pair production on the (poorly known 
H) extragalactic radio background; the mean free path is estimated to be ~ 1 — 5 Mpc at 10 11 GeV 0. 

^Further observations by Fly's Eye jl5| as well as its successor HiRES |l6| indicated a gradual increase 
in the depth of maximum with energy, suggesting a change in composition from heavy nuclei to nucleons. 
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In fact the possible existence of such high energy cosmic rays had already been mooted 
a decade earlier by particle physicists [[18] as a likely signature of Grand Unified Theo- 
ries (GUTs). In particular, Hill had proposed that relic topological defects such as 
magnetic monopoles created during GUT symmetry breaking in the early universe may 
have formed metastable bound states which would occasionally annihilate at the present 
epoch, thus releasing high energy quark and gluon jets which would hadronise to release 
cosmic rays of energy upto the GUT scale. Subsequently the development of superstring 
theory had led to independently motivated suggestions for the existence of massive relic 
particles associated with the 'hidden sector' of supersymmetry breaking; the lightest such 
state, having only gravitational interactions, would be metastable, thus a natural can- 



didate for the dark mater [20|. Indeed just before the Fly's Eye event was announced, 



detailed studies had been undertaken of the observational signatures of such particles |2T 



focussing on the possibility of detecting high energy neutrinos from their slow decays |2z 
(Moreover, the possibility that weakly interacting neutral particles such as neutrinos may 
acquire large interaction cross-sections through new physics had been discussed in the 
context of the detection of anomalous airshowers from the X-ray binary Cygnus X-3 ||23|| ). 

In recent years several other speculative ideas have been put forward for the possible 
origin of post-GZK UHECRs and some of them have been studied in sufficient detail as to 
allow clear experimental tests. The data set has also grown allowing some possibilities to 
be ruled out altogether and severe constraints to be put on others. This is an opportune 
time to assess the situation, particularly since the giant Pierre Auger Observatory ||24|| , 
under construction in Argentina, has recently seen "first light" and is expected to increase 
the world statistics of UHECRs by ten-fold within a few years. 



2 Observational Status 

The AGASA experiment has provided most of the data on UHECRs in recent years. 
In their 2000 review, Nagano and Watson [[l ]( had noted that the measurements were 
consistent at the level of ±15% in energy (or ±45% in flux) with older data from Haverah 
Park, Yakutsk and Fly's Eye. After making various necessary corrections, these authors 
had combined all the data to obtain the 'standard' differential energy spectrum: 

J ( E ) = C ( 63xl o9 G ev ) ~ 3 ' 20±0 ' 05 , for 4 x 10 8 GeV < E < 6.3 x 10 9 GeV (1) 

/ „ \ -2.75±0.20 n lr , 

= C ( 6,3x109 GeV ) > for 6 ' 3 X 10 GeV < E < 4 X 10 GeV > 

where C = (9.23 ±0.65) x 10 _33 m~ 2 s _1 sr _1 eV _1 . Although the spectral shape at higher 
energies was uncertain, there was no evidence of the GZK cutoff upto ~ 3 x 10 11 GeV. 



Preliminary data in 1999 from the HiRes air fluorescence detector |]2"5" |, the successor to 
Fly's Eye, had 7 events above 10 11 GeV, in conformity with this spectrum. 

At the International Cosmic Ray Conference in Hamburg in August 2001, AGASA 
reported 17 events above 10 11 GeV with a total exposure of ~ 6 x 10 16 m 2 sts (the data 
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set having been enlarged by accepting showers inclined by as much as 60° to the vertical), 
consistent with their earlier data p6[ . On the basis of more detailed comparison with 
simulations, the energy estimates were raised (by 20% at 10 11 GeV) relative to previously 
quoted results. As shown in Fig. |l], the flux at E > Eqzk is considerably above the ex- 
pectation for an uniformly distributed population of cosmologically distant sources (with 
a power-law injection spectrum chosen to match data at lower energies). 
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Figure 1: The AGASA spectrum compared with the expectation for the GZK-processed 
flux from cosmologically distant sources 



However HIRES, with a similar exposure in the mono mode, reported only 2 events 
above 10 11 GeV [^7[], instead of the ~ 20 events expected on the basis of the above stan- 
dard spectrum ([!]). In stereo mode, with about 20% of this exposure, only 1 event was 
seen. As seen in Fig. |2], there is a clear discrepancy with the AGASA data. The resolution 
of this may have to do with a possible underestimate of the efficiency of atmospheric fluo- 
rescence emission at longer wavelengths (more relevant for higher energies due to Rayleigh 
scattering) as suggested by a recent measurement PI . This may have caused HiRes to 
underestimate the energies (and overestimate the exposure) although the magnitude of 
such corrections cannot yet be precisely quantified [^9| . It is also the case that the begin- 
ing of the 'ankle' in the spectrum is a factor of ~ 2 lower in the fluorescence detector data 
than in the air shower detector data. The Auger Observatory employs both types of 
detectors so should be soon able to resolve these discrepancies and establish a consistent 
energy scale. However it should be emphasised that although the energy spectrum ([!]) 
may require revision, the existence of UHECRs above Eqzk is not in doubt! 
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Figure 2: New results from HiRes [27| compared with the 2-component spectrum su£ 
gested by the older Fly's Eye data . 



New results have also been presented by AGASA on the angular distribution of events 



on the sky |30fl . In accordance with earlier data, there is no evidence for any large-scale 
anisotropy for the 59 observed showers (having zenith angles < 45°) with energies above 
Eqzk- As seen in Fig. there is no observed correlation with the galactic plane (such 
as is seen at lower energies around ~ 10 9 GeV by both Fly's Eye and AGASA or 
with the super-galactic plane (as had been claimed earlier, particularly in the Haverah 
park data [0). However a number of 'clusters' — defined as a grouping of 2 or more 
events within (approximately the experimental angular resolution of) 2.5° — are seen; 
the chance probability that these result from an isotropic distribution is estimated by 
Monte Carlo to be 0.05% for the 5 observed doublets and 1.66% for the 1 observed triplet 
30fl . Assuming that all observed events above 4 x 10 10 GeV come from one particular 



class of source, the number of such sources is estimated assuming Poisson statistics to be 
220l?oo- It is also found that the cosmic rays contributing to the clusters have a hard 
energy spectrum with a differential spectral index of —1.8. If these are assumed to come 
from a separate source population then the total number of observed sources is limited 
to be less than 427 |30[]. However it has been pointed out that the data set used to make 
the initial claim for clustering ought not to be used in the actual analysis; moreover the 
directions of the most energetic events observed by previous experiments do not line up 



with any of the 6 'clusters' |29 |. Thus the indication for compact sources is exciting but 



still needs confirmation with increased statistics, as will be provided by Auger. 
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Figure 3: Arrival directions of UHECRs of energy E > 4 x 10 10 GeV (circles) and E > 
10 11 GeV (squares) [0. The area is weighted according to the AGASA exposure (the 
shaded area is unobservable). 



Finally, with regard to the composition of UHECRs, new data has muddled the picture 
suggested earlier by Fly's Eye and HiRes JO), [L6| of a change from heavy nucleus ( "iron" ) 
domination at 3 x 10 8 GeV to nucleon domination at 10 10 GeV. The muon content of 
the AGASA showers implies a mixed composition over ~ 10 9 — 10 10 GeV [[£J and a new 
analysis of Haverah Park data also indicates that nucleons make up no more than 34% 
of the flux in the range (0.3 — 2) x 10 9 GeV ||34|| . As seen in Fig.|], the conclusions are 
quite sensitive to the choice of the interaction model used [3j| hence the matter cannot 
be considered to be definitively settled as yet. However it would appear that a bound 
of < 50% can be set on the photon component at E > Eqzk from analysis of horizontal 
showers at Haverah Park, independently of the interaction model f34|. A similar limit 



is obtained by AGASA on the basis that photon-initiated showers tend to be muon- 
poor [ 36f| . This is potentially a powerful discriminant between models of the origin of 
UHECRs. Unfortunately the propagation of photons at these energies is uncertain since 
the dominant opacity (towards 77 pair production) is determined by the low frequency 
radio background which is poorly measured — the last experimental attempt was in 1970 
0. Subsequently mo del- dependent estimates have been made which suggest a mean free 
path at 10 11 GeV of a few Mpc |7J. Hence even if photons are absent in the UHECRs 
observed at Earth this does not exclude the possibility that they were emitted by the 
sources, especially if the mean free path turns out to be even smaller. 
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Figure 4: Variation of the depth of shower maximum with energy in different experiments, 
compared to expectations in various hadronic interaction models [[3fJ. 



3 Conventional explanations 

The lack of a GZK cutoff, which requires the sources to be relatively nearby, coupled with 
the isotropy, which implies a cosmologically distant population, pose severe problems for 
any astrophysical expalanation of UHECRs j37j. For example 7-ray bursts are isotrop- 
ically distributed and may well have the necessary capability to accelerate protons to 



such high energies ||38|| , but being cosmologically distant they cannot provide the events 
observed at E > E GZ k [0- Other possible sources such as active galactic nuclei may 
even be able to account for the post-GZK events if they are concentrated locally by a 
factor ^ 30 ||; however the local supercluster is overdense by a factor of only ~ 2. There 
are some nearby active galaxies but even if they can accelerate UHECR these would need 
to be deflected through large angles (and isotropised) by the intergalactic magnetic field 
(IGMF) |4(J; this requires the field strength to be ~ 100 times stronger than the upper 
limit of a few nG inferred from observational bounds on Faraday rotation in distant radio 
sources [[IT]. It has been proposed that the radio galaxy Cen A at 3.4 Mpc is the source 



of UHECRs however detailed calculations j|3| show that this requires an IGMF of 
~ 1 /iG to isotropise their directions and even so it would be difficult to account for the 
direction of the highest energy Fly's Eye event. By backtracing the 13 most energetic 
events in a magnetic field modelled as originating from a 'galactic wind', the radio galaxy 



M87 in the Virgo cluster has also been implicated as the source pHf . The convergence of 
the backtraced events turns out however to be due to the assumed magnetic field geometry 
and cannot therefore be considered evidence for an unique source [f4"5|. 
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4 Models involving new physics 



Given the problems of an astrophysical origin, it is natural to speculate that UHECRs are 
linked to physics beyond the Standard Model. Many novel suggestions have been made 
in this regard [JB|. I shall focus on only those which have been investigated in sufficient 



detail so as to provide falsifiability criteria. 
4.1 Z-bursts 

Since at least one species of relic neutrinos must have a mass of ~ 0.1 eV, it is an 
attractive possibility to suppose that these provide a target for ultra-high energy neutrinos 
from distant sources to annihilate on. This would create "Z-bursts" with an energy of 
m 2 z /2m l/ ~4x 10 12 (m !/ /l eV) -1 GeV, i.e. in just the right energy range to be a source for 
UHECRs |[47||. The energy spectrum of the nucleons and 7-rays resulting from Z decays 



is well measured at LEP so a detailed comparison can be made with the data f48] . In fact 
a good match to the energy spectrum is obtained for a relic neutrino mass of 0.26ig;? 4 eV; 
however the required ultra-high energy flux is rather large, taking into account that such 
light relic neutrinos cannot have a overdensity within the GZK distance of more than a 
factor of ~ 2 ||48|| . Although there are no direct experimental limits at such energies, any 
reasonable extrapolation to somewhat lower energies would violate the limits obtained 
recently by AMANDA. Moreover it is not clear how such high energy neutrinos would be 
created in the hypothetical cosmic sources — usually this would require the acceleration 
of even higher energy nucleons! (In principle, high energy neutrinos may be produced 
through the decays of super-massive relic particles |0| but such decays ought to also 



produce UHECRs directly as we discuss below, making this mechanism of less interest.) 
An alternative solution is to increase the local density of relic neutrinos by postulating that 



they are degenerate ||50|| ; however this is increasingly constrained by failure to observed 



the associated 'late ISW effect' in the angular power spectrum of the CMB, as discussed 



by Hannestad at this meeting |5T . 



4.2 Strongly interacting neutrinos 

Neutrinos with electroweak interactions cannot initiate the observed airshowers but, as 
has been speculated for some time |23[| , they may be become strongly interacting through 
new interactions. This possibility has been recently revived with the realisation that the 
string scale may in principle be as low as the TeV scale if there are new dimensions in 
Nature. However as discussed by Plumacher at this meeting, the increase in cross-section 
through (t-channel) exchange of Kaluza-Klein gravitons is inadequate to account for the 



characteristics of the observed events |p2f . A string-theoretic calculation in which the 
string amplitude is approximated as a sum over s-channel resonances (leptoquarks) also 
comes to a negative conclusion [^3| . 
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4.3 New neutral primaries 

The GZK cutoff would be higher if the primary particle is a new stable hadron which is 
heavier than a nucleon. For example it might be a bound state of a light gluino if that is 
the lightest supersymmetric particle |54j] or perhaps a H-dibaryon f55 [. The air showers 



that would be initiated by such a particle ("uhecron") have been studied in detail []5f| 
and comparison with data limits the primary mass to be ^ 50 GeV. The experimental 



constraints on new stable hadrons are already quite stringent |j7|] and it is likely that 
such a particle can soon be found or excluded. The interactions needed to produce such 
particles in distant sources must also produce neutrinos and 7-rays and the latter should 



be detectable, e.g. by Cerenkov telescopes such as VERITAS, MAGIC and HESS [46 



4.4 Violation of Lorentz invariance 

It has been noted that a small modification of the relation between momentum and energy 
in special relativity may undo the GZK cutoff f58|. In particular this may happen for a 
deformed dispersion relation c 2 p 2 = E 2 (l + E/Mp) + . . ., motivated by considerations 
of quantum gravity at the Planck scale f59fl . However this needs further examination in 
a quantum gravitational framework before definite conclusions can be drawn since, e.g., 



energy might be conserved only in a statistical sense |60 



4.5 Decaying supermassive dark matter 

As mentioned earlier, the possible existence of relic metastable massive particles whose 



decays can create high energy cosmic rays and neutrinos had been discussed |18, |2(J, 



before the famous Fly's Eye event [13| which focussed attention on the enigma of UHECRs. 
Subsequently this idea was revived f6l]] and it was further pointed out that such particles, 
being cold dark matter, would naturally have a overdensity by a factor of ~ 10 4 in the 
halo of our Galaxy |62|, |63[| . Hence if their slow decays generate UHECRs, all propagation 



effects will be unimportant (except possibly for photons) in determining the observed 
spectrum and cosmposition, and, moreover, there should be a detectable anisotropy |j4[ 
in the arrival directions, given our asymmetric position in the Galaxy. In order to account 
for the highest energy events with the observed rates, the particle mass should exceed 
mx ~ 10 12 GeV while its lifetime is determined to be tx — 3 x 10 9 Cx t , where to ~ 10 10 yr 
is the age of the universe and (x is the fraction of the halo dark matter density in the 
form of such particles []S5|| . This is in accordance with the theoretical expectation for 
'cryptons' [P0| . Moreover, as discussed by Kolb at this meeting, it is plausible that such 
particles ( "wimpzillas" ) were created with a cosmologically interesting abundance through 
the changing gravitational field at the end of inflation [^6| . 



The spectra of the decay products is essentially determined by the physics of QCD 
fragmentation ||65|| . Recently several groups have undertaken calculations which improve 
on previous efforts [J^. 63| in being both better connected to laboratory data (available 



S 



upto LEP energies) and incorporating the possible effects of new physics e.g. supersym- 
metry at higher energies [p7[ 



, ||, [TO], |7T], 0. As discussed by Toldra at this 
meeting, an effective approach is to evolve the fragmentation functions measured at the 
Z° peak (for neutrinos, inferred using HERWIG) upto the scale of the decaying particle 
mass by solution of the DGLAP evolution equations. Similar results are obtained in two 
independent calculations |)B, 7^] and, as seen in Fig. |5|, the evolved spectrum of nucleons 
is in good agreement with the 'flat' component ([[]) of cosmic rays at E > -Eqzk- On the 
negative side, the decay photons, which have a similar spectral shape are more abundant 
by a factor of ~ 2, so this model (as well as the related one involving annihilations of 
wimpzillas [Q) would seem to be ruled out by the Haverah Park bound |34| on the photon 
component of UHECRs. However, as commented earlier, it is quite possible that such high 
energy photons are significantly attenuated in their passage through the low frequency 
radio background; this would in turn generate a background of low energy 7-rays and it 
is necessary to check that this does not exceed observational limits (work in progress). 
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Figure 5: Fit to the UHECR spectrum beyond the 'ankle' with a decaying dark matter 
particle mass of 5 x 10 12 GeV (dashed line), including the effects of supersymmetry 
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Detailed calculations have also been made of the expected anisotropy, adopting differ- 
ent possible models of the dark matter halo (cusped, isothermal, triaxial and tilted) |75|] . 
We find that the amplitude of the anisotropy is controlled by the extent of the halo, while 
the phase is controlled by its shape. As seen in Fig. |6|, the amplitude of the first harmonic 
is ~ 0.5 for a cusped halo, but falls to ~ 0.3 for an isothermal halo, while the maximum 
is in the direction of the Galactic Centre, with deviations of up to 30° possible for triaxial 
and tilted haloes. Contrary to a recent claim |75| the halo of M31 is not bright enough 
to provide conclusive evidence for this hypothesis. However another key signature of a 
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dark matter origin may be some 'clustering' of events since the halo CDM distribution 
is expected to be clumpy rather than being smooth [|77j. Thus this model has several 
interesting observational signatures and can easily be falsified, e.g. if the UHECRs turn 
out to be heavy nuclei 178 . 




Figure 6: Contour plots (Hamer-Aitoff projections) of the UHECR sky for (a) cusped, 
(b) isothermal, (c) triaxial and (d) tilted models of the dark matter halo of our Galaxy. 
The effect of the halo of M31 is visible in the upper right of each plot |f75| . 



5 Conclusions 



After a long hiatus, high energy cosmic rays have again become very interesting for particle 
physicists looking for evidence of physics beyond the Standard Model. The source of the 
highest energy particles in Nature is an equally interesting enigma for astrophysicists. 
As Lemaitre first suggested, the origin of such particles may even be linked to the early 
universe, although not quite as he imagined! Presently the data are tantalising but not 
sufficient in either quantity or quality to distinguish definitively between proposed models. 
The good news that this will soon be remedied by the Auger array [Q] which should be 
complete by 2004 and is expected to detect ~ 2000 events above Eqzk within 5 years. 
Moreover the ambitious space-based experiment EUSO f7"9fl , scheduled for a 3- year flight 
on the International Space Station in 2007, will provide another substantial increase in 
collecting power. This is an exciting time for cosmic ray physics and we look forward to 
the surprises that Nature has in store for us. 
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